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SUMMARY 


An approximate solution for the potential flow through a cascade consisting of 
two rows of staggered circular rods is derived. The solution is then used to obtain 
the classic slotted-wall boundary-condition coefficient applicable to rod-wall wind 
tunnels. A comparison with the solution of Chen and Mears (as corrected by Barnwell) 
for flow through an unstaggered cascade is also made. 


INTRODUCTION 

Longitudinally slotted wind-tunnel test sections which have slat shapes with 
circular cross sections are generally referred to as rod-wall wind tunnels. Fig- 
ure 1 (a) shows the side view of the wind-tunnel flow expanding around the model and 
going through the wall into the plenum surrounding the test section. Figure 1(b) 
gives a cross-sectional view of the rod elements and defines the pertinent geometric 
parameters. 


Rod wall — y Plenum 

Tunnel 



(a) Side view of rod-wall wind tunnel. 



(b) Cross section A-A. End view of rod wall. 


Figure 1.- Diagram of rod-wall wind tunnel. 



In 1952, Sabol (ref. 1) published the ballistic range results of a study of the 
shock reflection properties of walls with various slat cross-sectional shapes 
(including circular) in both circular and rectangular test sections. His studies 
indicated that reflected shock waves of sharply reduced strength occurred when the 
walls had faired slots. In 1974, Binion and Anderson (ref. 2) corrpared acoustic data 
and aerodynamic data taken in both rod-wall and perforated-wall wind tunnels at 
transonic speeds. They found that the rod walls had lower noise levels than the 
perforated walls and that the aerodynamic data were as good as, or better than, the 
perforated-wall results. Gilliam (ref. 3) extended the transonic rod-wall acoustics 
study of reference 2 and found this type of wall to be very quiet in comparison with 
other transonic tunnels. His results also indicated that the subsonic fluctuating 
pressure coefficient approached that of a solid-wall wind tunnel with a turbulent 
wall boundary layer. Harvey et al. (ref. 4) examined the effect of rod gap width on 
noise attenuation as well as on wall boundary-layer removal and laminarization in 
supersonic tunnels. They were able to correlate the transition Reynolds number with 
gap width and to show that laminar rod-wall boundary layers introduced insignificant 
noise into the flow field. Recently, Harney et al* (ref. 5) published a paper citing 
their work on a transonic, adaptive rod-wall wind tunnel; it had solid sidewalls with 
upper and lower rod walls with both longitudinal and spanwise wall-jacking stations. 
This configuration allowed two-dimensional and some three-dimensional wall adaptation 
through changes in the wall openness ratio or changes in the wall shape or both to 
match free air streamlines about the model. 

The purpose of this report is to obtain an approximate, homogeneous rod-wall 
boundary condition. First, a solution for the inviscid, incompressible cross flow 
through a rod wall composed of alternately staggered rods is obtained by linear 
superposition of doublet singularities. After constructing this solution, it is then 
possible, using the method outlined by Barnwell (ref. 6), to obtain a classic homo- 
geneous s lo t te d-wa 1 1 -type boundary condition for the rod wall. This solution differs 
from that of Chen and Mears (ref. 7) in that their solution did not allow for off- 
axis movement of the rods. 


SYMBOLS 

a slot spacing (horizontal displacement between adjacent rods) 

c 1 doublet strength coefficient 

C 2 doublet-rod strength coefficient 

F complex potential functions 

h vertical displacement of adjacent rods 

Im( ) imaginary part of ( ) 

i imaginary number, pr 

K slotted-wall boundary -condition coefficient 

Re( ) real part of ( ) 
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rod radius 

horizontal cross-flow velocity 
vertical cross-flow velocity 

vertical cross-flow velocity far from the rods (i.e., the tunnel ambient 
component) 

complex velocity, U - iV 
cross-flow coordinates (see fig. 1 (b) ) 
coordinates in z'-plane (see fig. 2(a)) 
starting and ending locations of doublet rod 
coordinates for gap point 

vertical distance from rod center to isolated doublet 

stagnation point 

complex variable, x + iy 

intermediate complex variable, x' + iy' 

angles of inclination between doublet axis and x-axis 

slot width 

angle of inclination between doublet axis and y-axis 

transformation variable, i = 1, 2, 3, 4 

doublet strength 

strength of isolated doublet 

strength of doublet rod per unit length 

complex rod-wall potential, <t> + i4> 

nondimensional rapidly varying potential 

velocity potential 

velocity potential used to enforce far-field boundary condition 
rapidly varying portion of velocity potential 
stream function 
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FORMULATION AND ANALYSIS 


The problem presented herein is to suitably approximate the inviscid cross flow 
through a cascade of circular bodies as shown in figure 1(b). The potential is con- 
structed by adding combinations of elementary solutions of the Laplace equation which 
satisfy the boundary conditions at specific locations. It is assumed that no lift is 
generated; that is, the flow does not separate. Since the bodies are closed, only 
doublet singularities are needed. 

The complex potential of a doublet singularity at the origin of the z 1 -plane is 
given by 


F (z 1 ) 



( 1 ) 


If a row of doublets is placed at z* = 0, ±2a, ..., +2na, ... we have 


F (z 1 ) = -jie 


ia 


00 / \ 

, y ( i i \ 

z' ~ , Vz* + 2na + z* - 2na/ 

n=1 x 7 


Letting = z'/ a an ^ using egression (830) of reference 8 gives 


P(X 1 ) 
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Thus, 


F (z 1 ) 


ia 


e 7i 
“2a 



(2) 


If rows of doublets are placed at z* = z - iy d and z" = z + a - i(h - y d ) f then, 
from equation (2), the complex potential of the combined doublet rows in the z-plane 
is ejqpressed as 


F(z) 



(3) 


4 



where ^ represents the strength of each of the individual doublets and a = tc/ 2, 
which aligns the doublet axes with the oncoming flow through the wall. 

Now let equation (1 ) be the elemental potential of a continuous, constant- 
strength distribution of doublets, that is, n = |i 2 dx' . It can be seen by referring 
to figure 2(a) that 


F (z ' ) = -n 2 e 


io f a dx' IP p> 

L l 


dx' 


z' - x' 


ia , f z ' + x a\ ip , f z ' - x ' 

= »*2 e ln (7-Tl^j + »2 e ln ( — 


(4) 


where a = ( tc/2 ) + 0 and p = (n/2) - 0* If an infinite row of doublet-rod pairs 
are placed at z' , z 1 ± 2a, . then 
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Figure 2 



(b) Distributed singularities in z -plane. 

- Singularity distribution used for modeling the rod-wall 
boundary condition* 
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the complex potential function becomes 



after using expression (1016) of reference 8* Thus, 


ia , Sin 2T (Z + X a ip , Sln to (Z “ X b ) 

P(z- ) = |i e In — + In — — — " 

Sin 2T (Z + V Sln 2T (Z - X a 


Ihe doublet-rod pairs defined by equation (5) are next placed at 
z' = z - i(h/2) along with their images, which are placed at z - i(h/2) + a. Ihe 
combination of equation (5) plus its images is given by 
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The resulting system of singularities is as shown in figure 2(b). Equation (9) 
represents the complex potential of a cascade consisting of two staggered rows of 
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closed bodies whose shapes are determined by choosing , c 2 , h, y x , and 

Xjp* The complex velocity W is obtained by differentiating equation ( 9 ), which 
gives the following: 




sin 9 






( 10 ) 


The stream and potential functions may be obtained by taking the real and imaginary 
portions of equation (9). The velocity components may be obtained from the real and 
imaginary portions of equation (10). Thus, 


<j> = Re($) 

(1 la) 

cj> = Im(5>) 

(1 1b) 

U = Re (w) 

(1 1c) 

V = -Im(W) 

(1 Id) 
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which give the following: 


4) 
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v| 


= -y + c 1 A i 2 - =2 [sln 0 (s 21 • R 3l’ + COS 9 la 52 ' a 42 >] 
= X - o^,, + c 2 [sin 9 (A 22 - s 32 ) t cos 9 (A 41 - A &1 )] 
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z= iys 

z=iy * + °1 A 62 - °2 [sln 9 U 72 ‘ a 8 2> + COS 9 la 9 2 ‘ 


a i02>l - 


(1 2) 

( 13 ) 

( 14 ) 


( 15 ) 


where the A.^ coefficients are defined in the appendix and y g is the stagnation 
point. The problem is to now find or specify values of 9, c^, c 2 » y^, h/d, x a , 

and x b such that the desired rod shape is reasonably approximated. 


HOMOGENEOUS BOUNDARY CONDITION 
The classic form of the homogeneous boundary condition is 


d<t> a 2 <b 

■jT + aK -r — J- 
ox ox 5y 


0 


( 16 ) 


where K is the boundary -condition coefficient for a specified wall geometry. With 
Barnwell's method (ref. 6), K may be easily obtained once the potential is derived. 
Since solutions to Laplace's equation are correct to within an additive constant, 
equation (12) is written as 


<t> = -7 + 4> r + 4> 0 (17) 

where ^ is the rapidly varying portion of the potential in the vicinity of the 
wall and <t> Q is an additive constant which forces the boundary condition 

lim <t> = 0 (18 ) 

y-X© 


to hold. Thus, 


= lim ( <1> + y - <|> ) 
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( 19 ) 



giving 


♦o = ■ ? C 1 + 2c 2 (x b " X a ) COS Q ] 


( 20 ) 


From equation (20), it follows that 


= C l( A 12 + 7) + C 2 t sin 0 (A 


- A ) + cos 0 (A - A ) + J*(x K - x ) cos 0] 


21 31 


52 42' a b a ; 


( 21 ) 


Equation (21) may be written as 


d> = a<£ 
r r 


( 22 ) 


Evaluating 3?^ in the slot gives the desired result; that is, 

K = $ r (x = a/2; y = h/2) (23) 

or 

K = '^•[c 1 + 2c 2 (x b ~ x a * cos 0 -l (24) 

a 

As a check on the solution, equation (15) was compared with the Chen and Mears 
solution (ref. 7) as corrected by Barnwell (ref. 6). Chen and Mears used a doublet- 
rod model similar to what is presented here; however, their solution did not allow 
for the off-axis displacement of the alternate wall members. If c 1# x a > 0/ and 
h are set to zero, then equation (15) reduces to the corrected Chen and Mears 
solution. 


METHOD OF SOLUTION 

Ihe equations which must be solved are given by equations (13), (15), and (24). 
The unknowns are x a , x^, c-j , c 2 / 0/ and K, with parameters h, r, 

and a. Tb reduce the number of unknowns, x a is arbitrarily set to 0.01 for cases 
when h £ 0 and to 0 for cases when the rods are aligned. Ihe rod length is 

determined such that the zero streamline (eq. (1 3) ) is forced through a point in the 
gap with the coordinates 


y = h/2 (25a) 

g 


X 

g 



2 1/2 
- y ) 
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(25b) 
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It should be noted that the rod radius r sets the stagnation point y at x 
equal to zero. The doublet strength is determined such that equation (15) is 

satisfied at the upper stagnation point (0,y s ), and the doublet-rod strength c 2 is 
determined such that equation (15) is satisfied at the lower stagnation point 

<°,-y s ). 

In reference 6 it was shown that K was very dependent on the slot radius of 
curvature. Therefore, 0 is determined so that the radius of curvature in the gap 
is equal to r. The doublet position y d is set so as to minimize the root-mean- 
square (rms) deviation from a circular body of radius r. The boundary-condition 
coefficient K is calculated from equation (24). 

The solution procedure is iterative and, in many cases, the convergence is 
highly sensitive to the initial choice of variables. The gap point (x ,y ) where the 
radius of curvature is specified is not the point of minimum separation between the 
displaced cylinders. However, for small values of h, it is sufficiently accurate. 

If the minimum-separation point is used, instability develops in the solution pro- 
cedure and, for the larger wall openness values, undesirable body shapes are 
obtained. 


RESULTS 

Wall openness 6/a as a function of rod displacement h/a is given by the 
expression 


_6 

a 



(26) 


and is plotted for a rod radius of a/2 in figure 3. Body shapes for various values 
of wall openness and a fixed radius r of a/2 are shown in figure 4. These shapes 
were determined from equation (13) by solving for the zero streamline. The maximum 
rms error from a circular body for 0 < h/a < 0.90 was 0.0070, which occurred at 
h/a = 0.45. Typical rms values, however, were of the order of 0.0040. For those 
cases in which the rods were not offset, the openness was changed by varying r; the 
rms error decreased with increasing wall openness, going from 0.005 at 0.6-percent 
open to less than 0.001 at 30.0-percent open. 

In the present study, the rod radius and the radius of curvature matched in the 
gap between the doublet rods at the point (x , y ) as opposed to matching in the slot 
at the minimum. For openness ratios less than 10 percent, the difference in location 
between the minimum and the gap are insignificant. For openness ratios greater than 
10 percent, the difference can be substantial. Barnwell (ref. 6) showed large 
changes in the boundary-condition coefficient K with varying radius of curvature 
when the wall openness ratio was small. Large changes also occur with varying open- 
ness ratio and fixed curvature. During the present study it was found that these two 
parameters controlled the value of K for small openness ratios even more than the 
actual rod cross-sectional shape. For large openness ratios, the major effect is 
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10 




Figure 4.- Calculated body shapes for various values of h/a for rod 

radius of a/2. 


openness ratio. Figure 5 shews a comparison between the present case of displaced 
rods with a radius of a/2 and the Chen and Mears solution for undisplaced circular 
rods (ref. 7). For openness ratios of less than about 12 percent, the Chen and Mears 
solution and the present solution are equally applicable. If the openness exceeds 
about 12 percent, the effect of rod displacement begins to appear and the difference 
between the solutions grows to a factor of about 2. It should be emphasized, 
however, that the applicability of the homogeneous wall boundary condition becomes 
suspect at these large values of openness. 



Figure 5.- Rod-wall boundary condition for cylinder of radius a/2. 


CONCLUDING REMARKS 

An approximate, two-dimensional potential flow solution for the flow through a 
cascade of alternately displaced circular rods has been formulated. The solution was 
constructed by using discrete doublet and doublet-rod singularities. Free parameters 
in the solution were defined by minimizing the root-mean-square error deviation 
between the body streamline and a circle of given radius. Ihe radius of curvature in 
the gap between two adjacent bodies was also specified and matched. With a method 
described by Barnwell, the homogeneous slotted-wall wind-tunnel boundary-condition 
coefficient is derived. 
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For the case of undisplaced rods, the solution reduced to that of the classic 
Chen and Mears solution for the rod wall as corrected by Barnwell, For cases when 
alternate rods are displaced, the maximum root-mean-square error deviation from a 
circular cross section was determined as 0,0070, The boundary— condition coefficient 
is approximately equal to that given by the corrected Chen and Mears theory for wall 
openness ratio of less than about 12 percent. Values of wall openness ratio greater 
than this show deviations increasing to a factor of about 2. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 3, 1984 
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APPENDIX 


COEFFICIENTS OF STREAM, POTENTIAL, AND VELOCITY FUNCTIONS 


If we define the parameters 




P = ny d /a 


Q = Tlx /a 
a 


R = toc, / a 

D 


s ^ = Tix/a 


s = ti(x - x )/a 
^ a 


s = tx(x + x ) /a 
j a 


s = n;(x - x, )/a 
4 b 


s 5 = n;(x + x b )/a 


t 1 = it[y - (h/2) ]/a 


t 2 = n:[ (h/2) - Y d l/a ^ 


(A1 ) 


then the A^. coefficients used in equations (12) to (15) are given by the 
following: 


^ — sin 2s 1 - sinh t^ sin s 1 sinh t^ 

Hi = a" 2 2 

(sin s 1 cosh t 1 ) + (cos s 1 sinh t 1 + sinh t^) 


(A2) 


• 7 T sinh 2t + sinh t cos s cosh t 

v 4 1 2 1 1 

^12 — ~ a" 2 2 

(sin s 1 cosh ) + (cos s^ sinh + sinh t 


(A3) 
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APPENDIX 


2 2 

(-cos s cosh t + cos R) + (sin s sinh t ) 

A = 1 ln — 

212 2 2 
(cos cosh t ^ + cos R) + (sin s ^ sinh ) 


1 2 cos R sin s 1 sinh t * 

a 22 = tan " —r H- . 7T 

cos R - cos s 1 - smh 


2 2 

1 (-cos s^ cosh + cos Q) + (sin s ^ sinh ) 
31 “ (cos s 1 cosh t 1 + cos Q ) 2 + (sin s 1 sinh t ) 2 


_ 1 2 cos Q sin s 1 sinh t.. 

A 3 2 ' tan " T 2 — 2 “ 

cos Q - cos s ^ - sinh 


.2 , . , 2 , 

. sin s„ + smh t 
, 1 , 4 1 

A 41 = 2 ln — 2 —2— 

sm s c + sinh t 
o 1 


sinh 2t^ sin 2R 

^42 ^ an cos 2R cosh 2t^ - cos 2s 


2 2 
sin s + sinh t 

A 51 ’ 2 ln — 1 7T2 

sin s^ + sinh 


sinh 2 t 1 sin 2 Q 


A 52 tan cos cos h 2 t - cos 2 s 

1 1 


, _ x 2 1 - sinh — (y - ^-) sinh t„ 

a Vs 2/ 2 

\a/ r , . x -.2 


sinh 1 (y s - |) 


+ sinh t 
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APPENDIX 


u 2 sinh f (y s - 1) cos R 

* - 2 * - MD| (r. - |) 


2 sinh - 
it a 

a 2 n f 
cos Q - c 


. _ 2it sin 2R 

92 3 or, ^ 2it / 

cos 2R - cosh — y 

3 \ S 


(A13) 


(A14) 


(A15) 


sin 2 Q 


cos 2Q - cosh 


2 it / h\ 

i" ( y s - 2 ) 


(A16) 


Coefficients in equations (A1 ) to (All) are obtained from equations (11a) 
and (11b) and coefficients in equations (A12) to (A16) are obtained by evaluating 
equations (11c) and (lid) at x = 0 and y = y , the stagnation point on the 
cylinder centered at z = 0. 
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